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I Abstract. The bulk viscosities of two color-superconducting phases, the color-flavor locked (CFL) 

phase and the 2SC phase, are computed and compared to the result for unpaired quark matter. In 
D ■ the case of the CFL phase, processes involving kaons and the superfluid mode give the largest 

contribution to the bulk viscosity since all fermionic modes are gapped. In the case of the 2SC 
\^ • phase, ungapped fermionic modes are present and the process u + d ^ u + s provides the dominant 

' contribution. In both cases, the bulk viscosity can become larger than that of the unpaired phase 

for sufficiently large temperatures (T > I MeV for CFL, T > 0.1 MeV for 2SC). Bulk viscosity 
(as well as shear viscosity) is important for the damping of r-modes in compact stars and thus can 
potentially be used as an indirect signal for the presence or absence of color-superconducting quark 
matter. 
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Color superconductivity in compact stars. - Matter at sufficiently large densities 
and low temperatures is a color superconductor, which is a degenerate Fermi gas of 
if^ \ quarks with a condensate of Cooper pairs near the Fermi surface [1]. At asymptotically 

CN ■ large densities, where the quark masses are negligibly small compared to the quark 

chemical potential /i, three-flavor quark matter is in the color-flavor locked (CFL) 
^ ! state [0]. In this state quarks form Cooper pairs in a particularly symmetric fashion. 

I The gauge group of the strong interactions SU{3)c and the chiral symmetry group 

O ■ SU{3)i X SU{3)r are spontaneously broken down to the group of joint color-flavor 

rotations SU {3)c+l+r- All quarks participate in pairing, giving rise to energy gaps in 
K/i ; all quasifermion modes. Therefore, the lightest degrees of freedom in the CFL phase 

■ are Goldstone bosons: the broken chiral symmetry gives rise to an octet, with the 

lightest modes being the charged and neutral kaons (at not asymptotically large densities, 
chiral symmetry is an approximate symmetry and the Goldstone modes acquire a small 
mass). Moreover, the CFL phase is a superfluid, spontaneously breaking baryon number 
conservation symmetry U{1)b- Hence there is a "superfluid mode" which is exactly 
massless. Goldstone modes in the CFL phase can be described within effective theories 
[II, 0, 111], not unlike conventional chiral perturbation theory of the QCD vacuum. The 
results given below for the CFL bulk viscosity are computed within these effective 
theories. 

While rigorous QCD calculations from first principles can be done at asymptotically 
large densities, the situation is more complicated at moderate densities. Here we are 
interested in densities of matter inside a compact star. These densities can be as large 
as several times nuclear ground state density; however, even then the quark chemical 
potential is at most of the order of 500 MeV. Therefore, perturbative methods within 
QCD are not applicable and one has to rely on more phenomenological models. Fur- 



thermore, the mass of the strange quark is not negligible. This fact, together with 
the conditions of electric and color neutrality, imposes a "stress" on the CFL phase. By 
stress we mean here that the Fermi momenta of the quarks that pair are no longer equal. 
A sufficiently large stress (i.e., a sufficiently large mismatch in Fermi momenta com- 
pared to the pairing gap) will eventually disfavor the CFL phase; however, the density at 
which this happens is unknown. Unless the CFL phase is directly superseded by nuclear 
matter, less symmetric color- superconducting phases are expected to form the ground 
state. One of these phases is the so-called 2SC phase [HI]. In the 2SC phase, all strange 
quarks as well as up and down quarks of one color remain unpaired. This is one possible 
way for the system to deal with the heaviness of the strange quark. In phase diagrams 
based on Nambu-Jona-Lasinio models, the 2SC phase appears to be a candidate for the 
ground state at moderate densities laLZl]. Here we shall not discuss other interesting pos- 
sibilities such as Larkin-Ovchinnikov-Fulde-Ferrell (LOFF) states |8i, l9fl or spin-triplet 



pairing [|10l . llin which break translational and/or rotational invariance. For a calculation 
of the bulk viscosity in spin-triplet phases, see Ref. fl^. 

The only system in nature where deconfined quark matter may be expected, and 
hence the only "laboratory" for the study of color superconductivity, is the interior of 
compact stars. The temperature of compact stars drops below ~ 1 MeV within the first 
minute of the star's life. Hence, quark matter, if present in the core of the star, can be 
expected to be in a color- superconducting state, since the critical temperature for spin- 
singlet color superconductors, such as the 2SC and CFL phases, is estimated to be of the 
order of tens of MeV. Consequently, it is of great interest to compute thermodynamical 
and transport properties of different candidate color superconductors. Results of these 
calculations can be used to predict the behavior of the star with a quark matter core in the 
respective phase. Comparing these predictions with the actual astrophysical observations 
then may serve to rule out certain possibilities for the ground state of quark matter and 
thus eventually lead to a better understanding of the phases of QCD at high densities. 
One observable which is sensitive to the micsoscopic properties of the matter inside 
the star is for example the temperature as a function of time, i.e., the cooling curve 
of a star. Since the dominant cooling mechanism within the first million years in the 
evolution of a compact star is provided by energy loss through neutrinos, calculations 
of the neutrino emissivity (and the specific heat) have been performed for a variety of 
color superconductors, see for instance [Tsi 14, Ts]]. The remainder of these proceedings 



is devoted to another transport property, namely bulk viscosity, which is related to the 
rotation frequency of the star. 

Bulk viscosity and r-modes. - A compact star exhibits vibrational or rotational modes. 



typically in the kHz range [|16l] . The spectrum of these pulsations is very rich. They have 
been classified according to their restoring force. For instance, so-called r-modes are 
associated with the coriolis force. These modes are of particular interest because they 
are unstable with respect to emission of gravitational waves for all perfect fluid stars 



[|17LI18|] . For a nice pedagogical introduction on neutron star pulsations and instabilities, 
see Ref. If a neutron star was a perfect fluid, r-mode instabilities would lead to 
a drastic and quick spin-down, since the rotational energy would be transformed into 
the energy of the gravitational waves. The fact that fast rotating stars are observed leads 
to the conclusion that there must be a damping mechanism for these instabilities. This 



damping can be provided by the viscosity of the fluid. 

Both shear and bulk viscosity can provide damping. In Fig. [U we show that they are 
important in different temperature regimes in the example of unpaired quark matter. In 
the following, we discuss if and how color superconductivity may change this picture. 



mm. 



Calculations of bulk viscosity exist for conventional neutron star matter [|2( 

hyperonic matter [[23l. l24l] . and kaon condensed hadronic matter jlS] . For bulk viscosity 
in unpaired quark matter, see Refs. ll26l. l27ll28l[29L [30ll. 
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FIGURE 1. Damping effects of shear and bulk viscosity of unpaired quark matter in the rotation 
frequency - temperature diagram for a strange mass = 200 MeV. Shaded areas are forbidden due 
to r-mode instabilities. The damping effects are dominated by shear viscosity for temperatures lower than 
T ^ IQ^K, while they are dominated by bulk viscosity for larger temperatures, 10^ K < T < IO^'K. The 
dotted lines correspond to a different strange mass, M, = 100 MeV. The dashed line in the middle of the 
diagram corresponds to nuclear matter ms = millisecond pulsars, LMXB = low-mass X-ray binaries. The 
figure is taken (and slightly modified) from Ref. 131,1 . 



Definition of bulk viscosity. - The bulk viscosity ^ is defined as 



(1) 



It is proportional to the average dissipated power per volume {dW /dt)'m one oscillation 
period x = 2tz/o). In the present case, (local) volume oscillations are created by r- 
modes, whose typical frequencies are of the order of the rotation frequency of the star, 
®/ (27r) < Ims^^ (Another application of bulk viscosity is radial pulsations of neutron 
stars ll32U 33n.) We denote the volume oscillation by V{t) = Vq + 5Vqco% cot. In general, 
a superfluid has more than one bulk viscosity coefficient, since there are contributions 
related to the stresses in the superfluid flow relative to the normal one ll34l. [35l [36ll. We 



shall neglect these effects and compute the single bulk viscosity given in Eq. ([T]). (In the 
usual terminology of superfluid hydrodynamics our ^ is equivalent to the coefficient ^2-) 
Let us first consider a system where the lightest degrees of freedom are fermions, 
e.g., unpaired and 2SC quark matter. The volume oscillation forces the system out of 
chemical equilibrium with respect to the nonleptonic process 



and the semi-leptonic processes 



u + d ^ u + s (2) 

u + e ^ d + Ve (3a) 
u + e ^ s + Vg. (3b) 



In other words, a volume oscillation causes oscillations in the quantities = Us — lid, 
5/i2 = /irf — llu— lie, §113 = lis — 111, — lie (ouly two of which are independent). In the 
following, we ignore the semi-leptonic processes, which are in general slower than 
the nonleptonic one, and denote 8ll = The effect of the semi-leptonic processes 



has been discussed in Refs. pOL I37|] . The weak interaction provides the dominating 
processes since they "resonate" best with the external oscillation frequency. In principle, 
strong interactions also provide a contribution to the bulk viscosity, e.g., through thermal 
reequilibration. However, the speed of these processes compared to the time scale set by 
the oscillation renders this contribution negligible. 

The dissipated power is related to the oscillations in volume and chemical potential, 

BX{5n{t)5v{t)). (4) 



/dW\ 
YdT/ 



The coefficient X has to be determined from the microscopic calculation. It is defined 
by r^[5/i(?)] = A 5il{t), where rrf[5/i(?)] is the number of d quarks per volume and 
time produced in the reaction u + s ^ u + d. The linear dependence of on 5/i is a 
valid approximation for sufficiently small volume oscillations, 5Vo -C Vq. The coefficient 
B is related to the changes of the equilibrium particle densities with respect to the 
corresponding chemical potentials and depends on the phase. We have abbreviated 
dv = 5Vo/Vo cos cot. Finally, the relation between 5/i and 6v is given by the differential 
equation 

ddll ^ ddv „ , , _ 

where 7= CA is the characteristic inverse time scale of the microscopic process. The 
coefficient C depends, as the coefficient B above, only on the equilibrium densities and 
has to be determined for each phase separately. As pointed out in Ref. [[37|] . Eqs. dU 
and © show that the present system is completely analogous to an electric circuit with 
alternating voltage - playing the role of 5v{t) - and alternating current - playing the 
role of 6il{t). In this analogy, resistance R and capacitance C' can be identified with 
microscopic quantities of the quark system, l/B,C' ^ B/y. 

Solving the differential equation Q and inserting the result into Eq. (U) and then into 
([B yields 

B^ 7 



This expression for the bulk viscosity shows that, for a given ft), ^ has a maximum at 
7 = ft). Consequently, bulk viscosity is a resonance phenomenon with the largest bulk 
viscosity occurring when the time scale of the microscopic process matches the one of 
the external oscillation. 

In the CFL phase, the contributions to the bulk viscosity from the processes ^ and 
([3]) are exponentially suppressed and thus do not play any significant role for the bulk 
viscosity at any temperature ll3lll . The dominant contribution rather comes from the 
lightest Goldstone modes, the superfluid mode (p, originating from breaking baryon 
number conservation symmetry, and the neutral kaon K^, originating from the breaking 
of chiral symmetry. The contribution to ^ from the process 

9^9 + 9 (7) 

has been computed in Ref. While the (p is exactly massless, the mass of the 
can be deduced from the effective theory. The result, however, depends on the quark 
masses, the fermionic energy gap, and the quark chemical potential and is, strictly 
speaking, only valid at asymptotically large densities. Therefore, the (zero-temperature) 
kaon mass m^o (and the effective kaon chemical potential /i^o) are poorly known. In 
particular, it is not known whether the kaon mass and chemical potential allow for 
Bose-Einstein condensation. Condensation occurs for /i^^o > m^o, which seems likely 
from extrapolating the high density results to densities relevant for compact stars. If 
indeed the inequality /i^o > m^j-o is true, it can be expected that a kaon condensate 
is present for all relevant temperatures, since the critical temperature for condensation 



within the effective theory has been estimated to be of the order of tens of MeV 113911 . 
In the following, however, we shall assume that jlf^o < m^o, i.e., there are only thermal 
kaons in the system. The bulk viscosity in the presence of a condensate has not yet been 
computed. 

The contribution of the process 

K^^(p + (p (8) 

has been studied in Ref. [|40|l . This process reequilibrates the system with respect to the 
"out-of-equilibrium" quantity = /i^ — /i^, since K'^ carries quantum numbers sd. 

Rate of nonleptonic processes. - Before we present the results for the bulk viscosity 
for the 2SC and CFL phases, we briefly outline the calculation of the rate of the 
nonleptonic process (O. This calculation is done with a fixed nonzero yielding 
different rates for the two directions u-\-d ^u + s and u-\-s ^ u + d. Hence one obtains a 
net production rate of d quarks. In the 2SC phase, F^ is computed from the diagrams 
(and their reverse directions) shown in Fig. [21 The combinatorical factors in front of 
the diagrams are obtained upon counting color degrees of freedom: one can attach one 
of three colors to each of the two weak vertices, giving rise to 9 possibilities. In the 
2SC phase all blue quarks and all strange quarks are unpaired while all other modes are 
paired. Consequently, 4 of the 9 possibilities contain three gapped modes (red or green 
for both vertices), 2 contain two gapped modes (red or green for one, blue for the other 
vertex), 2 contain one gapped mode (blue for one, red or green for the other vertex), and 
one contains only unpaired modes (blue for both vertices). Therefore, at temperatures 





FIGURE 2. Contributions to the process u + d ^ u + s'm the 2SC phase. A gapped fermion is marked 
with the gap A at the respective line. We have omitted (small) contributions from anomalous propagators. 



much smaller than the fermionic energy gap, T ^ A, where the contributions of gapped 



quarks are exponentially suppressed, is to a very good approximation given by Bill 



rfc = lr™p forr<A, (9) 



since only the one reaction containing only unpaired modes contributes. The rate 



lUnp 



was computed in 114 111 . For larger temperatures, the contributions from the gapped modes 
cannot be neglected. They become important for temperatures 0.3 Tc <T<Tc. Each 
diagram yields a contribution which, for one direction of the process, schematically 



reads (for the exact expression and its derivation see Ref. [|37l1 ) 



X f{eiei)f{e2£2)f{-e3£3)f{-e4£4) ■ (10) 

Here, ^ is a complicated function of the momenta p, and the signs e, = ±1, / < 4; £, are 
the quasiparticle energies, and / the Fermi distribution functions. The sum over the signs 
ei is t ypi cal for collision integrals for any superfluid or superconductor (see for instance 
Ref. 1I42I1 in the context of superfluid ^He): the process u + d ^ u + s does not only 
get contributions from two quasiparticles u, d, going into two quasiparticles u, s. There 
are also contributions from three particles coalescing into one particle (and one particle 
decaying into three) and from annihilating (and creating) four particles. This is possible 
because quasiparticles can be created from or deposited into the condensate. Formally, 
this is obvious from the Fermi distribution functions where the factor /(e) accounts for 
an ingoing particle and /(— e) = 1 — /(e) for an outgoing particle. Due to the sum over 
{ei}, all combinations of the /'s can appear, e.g., /(ei)/(e2)/(e3) [1 — /(e4)]. Here we 
shall not present the explicit results for the rate itself but rather directly skip to the result 
for the bulk viscosity, which is a function of this rate. 

Results for 2SC and CFL quark matter - Inserting the result for the net production 
rate of d quarks into ^ yields the bulk viscosity for the 2SC phase as a function of 
temperature. The result for a typical oscillation frequency t = 2;r/(U = 1 ms is shown 
in Fig. |3l A critical temperature of Tc = SOMeV is assumed. For low temperatures, the 
time scale of the nonleptonic process is much smaller than the oscillation frequency 
Y <^ 0), implying C °^ 7- Consequently, from Eq. ^ we conclude ^2SC = Cunp/9. For 
large temperatures, however, we have 7 ^ ft) and thus ^ oc l/y. Consequently, the su- 
perconducting phase, which has the slower rate, has the larger bulk viscosity. The result 



becomes nontrivial (= ^2SC not being simply a numerical factor times ^unp) between 
these two regimes, i.e., 7~ ft), and for temperatures close to the critical temperatures, 
not shown in this figure (see Ref. 113711 ). 

The bulk viscosity from the process ([8]) in the CFL phase is obtained along the same 
lines. One inserts the result for the rate of the kaonic process dH) into Eq. Q. The result 
is shown in Fig. [3] for different values of 5m = m^o — jl^o (remember that we assume 
that there is no kaon condensate, i.e., dm > 0). We have also plotted the contribution 
from the process (|7]), which, in the limit ft) — > is given by [|38ll 



CFL 



0.011-^. 
T 



(11) 



We observe that at T ~ ( 1 — 10) MeV the bulk viscosity of CFL matter is comparable 
to that of unpaired quark matter. For T < I MeV, ^ is strongly suppressed. Depending 
on the poorly known value for dm the pure (p contribution (fTTI) may dominate over the 
contribution from the ^ (p + (p reaction at low enough temperatures. However, for 
r < 0. 1 MeV the (p mean free path is on the order of the size of the star, i.e., the system 
is in the coUisionless rather than in the hydrodynamic regime, and the result ceases to be 
meaningful. 
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FIGURE 3. Bulk viscosities as functions of temperature for an oscillation period i ~2n/(£)^ 1 ms. 
CFL phase: contribution from the process K'-^ ^ (p + (p for different values of 5m = niffi — jif^o and 
contribution from (p ^ (p + (p, see Eq. (fTTT i. 2SC phase and unpaired quark matter: contribution from 
the process u + d ^ u + s. 
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